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ABSTRACT: Real-time quantitative PCR (qRT-PCR) is one of the important methods for investigating the changes in mRNA 
expression levels in cells and tissues. Selection of the proper reference genes is very important when calibrating the results of real-time 
quantitative PCR. Studies on the selection of reference genes in goat tissues are limited, despite the economic importance of their meat 
and dairy products. We used real-time quantitative PCR to detect the expression levels of eight reference gene candidates (18S, TBP, 
HMBS, YWHAZ, ACTB, HPRT1, GAPDH and EEF1A2) in ten tissues types sourced from Boer goats. The optimal reference gene 
combination was selected according to the results determined by geNorm, NormFinder and Bestkeeper software packages. The analyses 
showed that tissue is an important variability factor in genes expression stability. When all tissues were considered, 18S, TBP and HMBS 
is the optimal reference combination for calibrating quantitative PCR analysis of gene expression from goat tissues. Dividing data set by 
tissues, ACTB was the most stable in stomach, small intestine and ovary, 18S in heart and spleen, HMBS in uterus and lung, TBP in liver, 
HPRT1 in kidney and GAPDH in muscle. Overall, this study provided valuable information about the goat reference genes that can be 
used in order to perform a proper normalisation when relative quantification by qRT-PCR studies is undertaken. (Key Words: Reference 
Gene, Expression Stability, Tissue, Goat) 



INTRODUCTION 

Goat is an economically important animal in many 
communities worldwide, providing meat and milk products 
with abundant nutritional value (Dubeuf et al., 2004). 
Previous studies have shown that breeding, growth and 
other major economical traits of goats are highly regulated 
by functional genes. Studies involving temporal and spatial 
expression specificity of these genes in different goat 
tissues are helpful when screening the functional genes 
regulating important economic traits of goats. They are also 
vital in marker-assisted selection (MAS) studies and in 
breeding program improvements (Moioli et al., 2007; 
Domby et al., 2010; Fu et al., 2011; Zi et al., 2012). 

Real-time quantitative PCR is a quantitative analytical 
method used to determine mRNA expression at the 
transcriptional level. It is rapid, highly reproducible and 
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extremely sensitive, and is therefore widely used in the 
quantitative studies on gene expression (Dabek et al., 2010; 
Novinscak et al., 2011; Veys et al., 2012). However, the 
reliability of the results from qRT-PCR is affected by 
experimental conditions, yield and quality of the sample 
RNA, as well as the efficiency of reverse transcription 
(Fleige and Pfaffl, 2006). Therefore, specific reference 
genes are normally selected for calibration and 
normalization of target genes in the studies in order to 
reduce the deviation of the results (Mehta et al., 2010; 
Wang and Xu, 2010; Pierzchala et al., 2011). Studies on the 
screening of reference genes used in qRT-PCR have been 
reported in tissue sourced from humans, mice, pigs, cattle 
and other mammals (Garcia-Crespo et al., 2005; Spalenza et 
al., 2011; Chechi et al., 2012; Zhang et al., 2012), but 
studies in goat tissues are limited. The expression levels of 
frequently used reference genes may differ significantly 
among different species, cell types, tissues, different stages 
of cell proliferations and organ development, in vitro 
culture and other kinds of experimental conditions (Chechi 
et al., 2012; Chi et al., 2012). Therefore, it is of great 
importance to select suitable reference genes to carry out 
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calibration and normalization for the expression of target 
genes in the specific experimental conditions. 

We selected eight candidates of reference genes (18S, 
TBP, HMBS, YWHAZ, ACTB, HPRT1, GAPDH, and 
EEF1A2), and the differences in their expression levels in 
ten different tissues types (heart, liver, spleen, lung, kidney, 
stomach, uterus, ovaries, small intestines and muscles) of 
Boer goat were compared. geNorm, NormFinder and 
Bestkeeper analytical software packages were used to 
determine the gene combination showing the most stable 
expression. Identifying the best combination of reference 
genes will provide information to select the most suitable 
reference genes for calibrating qRT-PCR analysis of 
mRNAs from goats. 

MATERIALS AND METHODS 

Tissues 

Ten different goat tissues (heart, liver, spleen, lung, 
kidney, stomach, small intestine, uterus, ovary and skeletal 
muscle) were collected from three 5 -month-old female Boer 
goats. All the samples were immediately snap-frozen in 
liquid nitrogen and stored at -80°C and would be used for 
evaluating the stability of each candidate genes. 

Total RNA isolation and cDNA synthesis 

Total RNA was extracted using Trizol (Invitrogen, 
Gaithersburg, MD) according to the manufacturer's 
instructions. RNA concentration was determined by 
spectrophotometry at 260 nm/280 nm absorbance ratio with 
expected values between 1.8-2.00 using a Nanodrop ND- 
1000 spectrophotometer (Nanodrop Technologies, 
Wilmington, DE, USA). RNA integrity was confirmed by 
gel electrophoresis using 1% agarose with ethidium 
bromide. RNA samples were run on a 1% agarose gel at 
100 V for 30 min (Supplementary Figure 1). The isolated 
RNA was immediately frozen at -80°C. cDNA was 
synthesized from 1 |ig of total RNA by using a PrimerScript 
RT reagent kit (TaKaRa, Osaka, Japan) with an oligo-dT 
primer. The mixture was incubated at 37°C for 15 min and 
85°C for five seconds to terminate reaction in an ABI 9700 



DNA Thermal Cycler (Applied Biosystems, Foster City, CA, 
USA). 

Reference gene selection and primer designing 

Eight genes were selected to be evaluated as reference 
genes according to the literature: 18S, TBP, HMBS, YWHAZ, 
ACTB, HPRT1, GAPDH, and EEF1A2 (Garcia-Crespo et al., 
2005; Spalenza et al., 2011; Chechi et al., 2012; Zhang et al., 
2012). Selection was done considering their stability values 
from published studies. The candidate genes and primers 
are listed in Table 1. The primers were designed using 
Primer Express 2.0 software (Yoon et al., 2002) and 
synthesized by Invitrogen Shanghai (Shanghai, China). 

Real-time quantitative PCR 

Real-time quantitative PCR was performed using three 
biological replicates and technical triplicates of each cDNA 
sample, according to the manufacturer's protocol. cDNA 
samples were made into seven dilution series (1:10), and 
each dilution was amplified in triplicate PCR amplifications 
and plotted as mean values to generate a standard curve. 
PCR was performed using the same volume of cDNA 
sample for each gene to record the cycle threshold (Ct) 
value. PCR experiments were performed with an ABI Prism 
7900 HT sequence-detection system (Applied Biosystems), 
using SYBR Green PCR Master Mix according to the 
manufacturer's protocol. A 50 uL reaction sample included 
5 |iL SYBR Green Real-time PCR Master Mix, 0.2 |iL (10 
|iM) primers, 1 |iL sample, 0.2 |iL ROX reference dye, and 
3.4 |iL nuclease-free water. The following cycling 
conditions were used: 95°C for 30 s and 40 cycles of 95°C 
for five seconds, and 60°C for 30 s. Fluorescence signal 
was obtained at a temperature of 72°C extension. 

Data analysis 

All data are representative of experiments performed at 
least three times in duplicate. Data are meanistandard 
deviation (SD). The Ct value of each gene after qRT-PCR 
was input into an Excel datasheet (Microsoft Excel 2003) 
before importing the data into geNorm (Vandesompele et al., 
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Supplementary Figure 1. Result of the gel electrophoresis of total RNA isolated from goat tissue samples. 1: heart, 2: liver, 3: spleen, 4: 
lung, 5: kidney, 6: stomach, 7: small intestine, 8: uterus, 9; ovary, 10: muscle. RNA samples were run on a 1% agarose gel at 100 V for 
30 min. All the samples show high-quality RNA because there is a clear appearance of the 28S, 18S and 5S rRNA bands. The results 
indicated that total RNA was isolated successfully and with no degradation and could be used for further research. 
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Table 1. Selected candidate reference genes used in the qRT-PCR assay indicating name, GenBank accession number, annealing 
temperature (T.J, PCR efficiency, slope, regression coefficient and primers used for the expression study 



Gene 



Full gene name 



GenBank 
accession 
number 



PCR 



Regression 



(°C) 



efficienc Slope coefficient Forward primer (5 '-3') 



Reverse primer (5 '-3') 



y (%) 



ACTB 
CAPDH 



TBP 
HPRT1 



YWHAZ 



HMBS 



I8S 

EEF1A2 



Beta-actin DQ845171 
Glyceraldehyde-3-phosphate AF017079 

dehydrogenase 
TATA box-binding protein DQ845178 
Hypoxanthine DQ845175 

phosphoribosyltransferase 1 
Tyrosine 3- 

monooxygenase/tryptophan 

5-Monooxygenase activation 

protein zeta polypeptide 
Hydroxymethylbilane synthase 

mRNA 

18 S ribosomal RNA DQ22245 3 . 1 

Eukaryotic translation BC108110.1 
elongation factor 1 alpha 2 



62 99.66 
62 96.84 



62 94.92 
60 102.21 



DQ845179 60 98.84 



DQ845174 60 97.63 



60 95.68 
60 101.78 



-3.33 
-3.40 



-3.45 
-3.27 



-3.35 



-3.39 



-3.43 
-3.29 



0.998 TGCCCTGAGGCTCTCTTCCA 
0.998 GCAAGTTCCACGGCACAGTC 



0.993 AAAGACCATTGCACTTCGTGC 
0.998 ATTATGGACAGGACCGAACGG 



TGCGGATGTCGACGTCACA 
CCCACTTGATGTTGGCAGGA 



ATATCAGTGCGGTGGTGCG 
CCAACAGGTCGGCAAAGAACT 



0.999 CCAACGCTTCACAAGCAGAGA TGCTTGCTGTGACTGGTCCA 



0.998 AGCATGCCTTGGAGAGGAATG CGCTTGCAGACAGCTCCAAT 



0.993 CGGCTACCACATCCAAGGAA 
0.997 CGACTGGCCACCTCATCTACA 



GCTGGAATTACCGCGGCT 
CCAGGCATACTTGAAGGAGCC 



2002), NormFinder (Pfaffl et al., 2004) and Bestkeeper 
(Andersen et al., 2004) software packages. 

geNorm was used to evaluate the stability of reference 
genes and for screening the optimal combination by using 
two parameters M (average expression stability) and V 
(pairwise variation). geNorm generates an M value for each 
gene that is arbitrarily suggested to be lower than 1.5 (with 
a lower value indicating increased gene stability across 
samples), and a pairwise stability measure (V number) to 
determine the benefit of adding extra reference genes for 
the normalization process (again with a lower value 
indicating greater stability of the normalization factor). 
Lower M values for the references indicated higher stability. 
Paired variation analysis was carried out using geNorm 
after the optimal combination of two reference genes was 
obtained. If the paired variation coefficient was lower than 
0.15, the combination of these two genes was sufficient to 
be used as the criteria for calibrating the results from 
quantitative analysis, and the next gene was not required. 
By contrast, if the paired variation coefficient was higher 
than 0.15, the next genes should be added to the analysis 
(Mehta et al., 2010). 



Similarly, NormFinder generates a stability measure of 
which a lower value indicates increased stability in gene 
expression and group samples. This allows a direct 
estimation of expression variation, ranking genes according 
to the similarity of their expression profiles by using a 
model-based approach. Bestkeeper generates a pairwise 
correlation coefficient between each gene and the 
Bestkeeper index (the geometric mean of the threshold 
cycle values of all the reference genes grouped together) 
(Perez et al., 2008). 

RESULTS 

Amplification efficiency and specificity of qRT-PCR 

As shown in Table 1, the regression coefficients of the 
standard curves for the eight candidates reference genes 
were all higher than 0.99. The amplification efficiencies 
were all higher than 94% and showed good linear 
correlation, indicating that the specificity of the primers was 
relatively good and that the results of the quantitative 
analysis were accurate and reliable. Mean Ct values and 
standard deviation for the genes in each tissue are shown in 
Table 2. The obtained Ct and amplification efficiency 



Table 2. Mean cycle threshold (Ct) values and standard deviation for candidate reference genes in examined tissues 



Tissues 



ACTB 



GAPDH 



TBP 



HMBS 



HPRT1 



EEF1A2 



YWHAZ 



■ Mean±SD ■ 



18S 



Heart 


20.02±0.15 


18.7110.41 


28.5510.13 


25.5411.15 


24.6210.27 


17.8610.50 


20.8711.23 


10.3810.36 


Liver 


19.72±0.06 


20.48±0.43 


28.6510.90 


25.8510.82 


21.2412.13 


22.3110.21 


21.9010.89 


11.3611.50 


Spleen 


17.40±0.17 


21.3810.43 


28.2910.39 


27.3110.69 


23.7910.24 


21.8410.41 


20.6610.95 


10.7410.40 


Lung 


17.91±0.48 


22.1110.06 


28.0710.38 


26.8010.39 


23.6410.45 


21.6710.25 


20.4810.63 


10.9310.91 


Kidney 


18.86+0.45 


19.7710.37 


28.2010.25 


24.9210.16 


24.6310.07 


21.9810.12 


21.0010.31 


9.7710.39 


Stomach 


17.94±0.17 


21.1910.53 


28.1310.29 


25.4310.73 


22.3910.91 


21.4510.31 


18.8610.66 


10.2610.31 


Muscle 


20.56±0.68 


17.1910.69 


29.5110.14 


26.5811.04 


24.7210.69 


18.4710.73 


22.2611.22 


10.6710.38 


Small intestine 


18.03±1.01 


21.0010.91 


28.7910.82 


25.6911.40 


24.7210.62 


21.6711.06 


20.8211.44 


9.9910.18 


Uterus 


18.34±0.04 


20.9610.18 


27.7710.17 


26.21+0.16 


24.0610.22 


21.2010.11 


20.4810.75 


10.3910.16 


Ovary 


18.11+0.13 


19.9410.31 


27.2310.32 


26.1910.50 


23.5510.67 


20.6310.20 


19.9511.02 


9.8111.05 
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values were used in the geNorm, NormFinder and 
Bestkeeper algorithms to calculate expression stability of 
the investigated reference genes for each tissue. 

Stability of reference gene expression 

Firstly, the stability was evaluated taking into account 
the entire data (all tissues and goats) with geNorm, 
NormFinder and Bestkeeper algorithms. Table 3 lists the 
stability of reference gene expression obtained from the 
geNorm analysis (from high to low): 18S, TBP, YWHAZ, 
HMBS, ACTB, HPRT1, GAPDH, and EEF1A2 (Figure 1). 
The stability of the reference gene expression obtained by 
NormFinder analysis is as follows (from high to low): 18S, 
TBP, HMBS, YWHAZ, ACTB, HPRT1, GAPDH, and 
EEF1A2. The stability of the reference gene expression 
obtained using Bestkeeper analysis is (from high to low): 
18S, TBP, HMBS, YWHAZ, HPRT1, ACTB, GAPDH, and 
EEF1A2. 

From the analytical results of the three software 
packages, it can be deduced that 18S and TBP displayed the 
best gene expression stability, while GAPDH and EEF1A2 
are the least stable. The gene expression stability of HMBS, 
YWHAZ, ACTB, and HPRT1 is moderate. 

The paired variation coefficient of geNorm was higher 

1.6 
| 1.4 
« 1.2 
B 1.0 
> 0.8 
£ 0.6 
I 0.4 
£ 0.2 

0.0 



Table 3. Expression stability values of the candidate reference 
genes calculated by geNorm, NormFinder and Bestkeeper 
algorithms (ranking in parentheses) 



Gene 


Stability value 


Stability value 


Stability value 


geNorm 


NormFinder 


Bestkeeper 


18S 


0.634 (1) 


0.156(1) 


0.396(1) 


TBP 


0.634 (2) 


0.467 (2) 


0.445 (2) 


YWHAZ 


0.722 (3) 


0.813 (4) 


0.642 (4) 


HMBS 


0.799 (4) 


0.584 (3) 


0.566 (3) 


ACTB 


0.887 (5) 


1.203 (5) 


0.881 (6) 


HPRT1 


1.048 (6) 


1.297 (6) 


0.825 (5) 


GAPDH 


1.297 (7) 


1.550 (7) 


1.096 (7) 


EEF1A2 


1.425 (8) 


1.585 (8) 


1.153 (8) 



than 0.15 (V 2/3 = 0.234) when 18S and TBP are used in 
combination as the reference genes, indicating that the 
reference gene combination with 18S and TBP was not 
sufficient to be used as the calibration criteria for the 
quantitative analysis for gene expression, and the next 
candidate for reference gene was required (Figure 2). When 
YWHAZ, HMBS, and ACTB (ranked in third, fourth and fifth 
places in expression stability) were successively added, the 
variation coefficient showed 0.197 (V 3 / 4 ), 0.183 (V 4/5 ) and 
0.210 (V 5/6 ). As 0.183 (V 4/5 ) was the smallest variation 



EEF1A2 GAPDH HPRT1 



ACTB 



HMBS YWHAZ 18SITBP 



*- Least stable genes Most stable genes — 

Figure 1. pGene expression stability of candidate reference genes. Gene expression stability of candidate reference genes in goat tissues 
analyzed by the geNorm program which proceeds to the stepwise exclusion of the genes whose relative expression levels are more 
variable among tissue samples. Threshold for eliminating a gene as unstable was M >1.5. Lower values of M correspond to the most 
stable genes, hence the most suitable for normalization. 



0.300 
0.250 

3 °' 20 ° 
1 0.150 

> 0.100 

0.050 

0.000 



0.234 



0.255 





0.206 




V2/3 



V3/4 V4/5 V5/6 V6/7 
Pairwise Variations 



V7/8 



Figure 2. Evaluation of the optimum number of reference genes according to the geNorm software. The magnitude of the change in the 
normalization factor after the inclusion of an additional reference gene reflects the improvement obtained. V ,*/;+/ represent the models 
being compared: those with ( and i+1 reference genes. 
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coefficient, so HMBS should be added and 18S, TBP, and 
HMBS should be the optimal combination of reference 
genes. 

The stability of the genes was also evaluated for each of 
the ten tissues analysed. As expected, the stability of the 
genes varies among tissues (Table 4). ACTB was the most 
stable gene in stomach, small intestine and ovary, 18S in 
heart and spleen, HMBS in uterus and lung, TBP in liver, 
HPRT1 in kidney and GAPDH in muscle, evidencing the 
specificity of each tissue developing characteristic 
biological functions and specific metabolic pathways. These 
results showed that tissue is an important variability factor 
that affects the stability genes expression. Thus, it is 
evidenced the necessity of using reference genes according 
the tissue analysed. 

DISCUSSION 

The techniques for quantitative analysis of gene 
expression at mRNA level include RNA blotting, gene 
microarray, and qRT-PCR. The final results from the 
expression of target genes should be calibrated using 
reference genes (Kidd et al., 2007; Kortner et al., 2011; 
Nesvadbova and Knoll, 2011; Lardizabal et al., 2012; 
Rivera- Vega et al., 2012). Ideal reference genes should be 
stably expressed in the tissues, cells and under the 
experimental conditions used in the experimental system. 
However, the so-called constant expression of any reference 
gene only refers to the specific conditions under certain 
environments, and the expression levels in different cell 
types, tissues and different experimental systems can be 
significant (Radonic et al., 2004; Bustin et al., 2009; Lanoix 
et al., 2012). Therefore, it is particularly important to use 
reference genes suitable for the specific experimental 
conditions and species under investigation. 

Studies on reference genes in mammals have previously 
been reported. Nygard et al. (2007) screened for reference 
genes in different tissues by using SYBR Green qRT-PCR 
technique. The results showed that the expression of ACTB, 



RPL4, TPB, and HPRT1 was the most stable; ACTB and 
RPL4 can be used as the reference genes for analyzing gene 
expression of high abundance genes, while TPB and HPRT1 
can be used as the reference genes for analyzing gene 
expression of low abundance genes (Nygard et al., 2007). 
Lisowski et al. (2008) screened reference genes in different 
tissues sourced from cattle. The results showed that the 
stability of the expression of reference genes was not the 
same in different tissues. The expression of ACTB and TBP 
was the most stable in liver, the expression of GAPDH and 
YWHAZ was the most stable in kidney tissue, the expression 
of GAPDH and SDHA was the most stable in pituitary 
tissue, and the expression of TBP and HPRTI was the most 
stable in thyroid (Lisowski et al., 2008). Crespo et al. 
(2005) found that different combinations from six reference 
genes (ACTB, YWHAZ, RPL19, GAPDH, G6PDH, and 
SDHA) were applicable for the calibration of the results 
from qRT-PCR analysis in different sheep tissues. For 
example, GAPDH, SDHA, and ACTB were suitable for 
cerebrum, G6PDH, ACTB, and YWHAZ were suitable for 
cerebellum, while SDHA, RPL19, and GAPDH were 
suitable for mesenteric lymph node (Garcia-Crespo et al., 
2005). As the optimal reference genes for different species 
and different tissues are not necessarily the same, it is 
essential to determine the optimal reference genes for 
specific species used in the experimental system. 

The present study was carried out to examine the 
expression levels of eight candidates reference genes (18S, 
TBP, HMBS, YWHAZ, ACTB, HPRTI, GAPDH, and 
EEF1A2) in heart, liver, spleen, lung, kidney, stomach, 
uterus, ovaries, small intestines, and muscles tissues 
sourced from Boer goat. Real-time quantitative PCR 
technique and three software packages (geNorm, 
NormFinder and Bestkeeper) were used to analyze the 
results. Though the algorithms utilized by the three software 
were slightly different, they all sorted the candidate 
reference genes by their level of expression stability. 

Even though the stability calculations of the eight 
reference genes obtained from geNorm, NormFinder and 



Table 4. Stability values (M-values) for each reference gene calculated by using geNorm algorithm 



Tissues 


ACTB 


GAPDH 


TBP 


HMBS 


HPRTI 


EEF1A2 


YWHAZ 


18S 


Heart 


0.390(5) 


0.348(4) 


0.260(3) 


0.853(8) 


0.433(6) 


0.173(2) 


0.671(7) 


0.171(1) 


Liver 


0.532(5) 


0.439(3) 


0.357(1) 


0.359(2) 


1.299(8) 


0.477(4) 


0.722(6) 


1.024(7) 


Spleen 


0.227(3) 


0.294(4) 


0.206(2) 


0.394(5) 


0.471(6) 


0.536(7) 


0.610(8) 


0.203(1) 


Lung 


0.413(6) 


0.354(4) 


0.311(3) 


0.242(1) 


0.245(2) 


0.387(5) 


0.623(8) 


0.518(7) 


Kidney 


0.261(7) 


0.239(6) 


0.176(4) 


0.148(3) 


0.051(1) 


0.055(2) 


0.324(8) 


0.221(5) 


Stomach 


0.148(1) 


0.407(5) 


0.301(3) 


0.512(6) 


0.640(8) 


0.152(2) 


0.577(7) 


0.373(4) 


Muscle 


0.714(8) 


0.147(1) 


0.620(6) 


0.381(3) 


0.556(5) 


0.149(2) 


0.430(4) 


0.664(7) 


Small intestine 


0.107(1) 


0.110(2) 


0.231(3) 


0.554(7) 


0.335(4) 


0.411(5) 


0.486(6) 


0.633(8) 


Uterus 


0.123(3) 


0.188(5) 


0.218(7) 


0.058(1) 


0.202(6) 


0.062(2) 


0.359(8) 


0.181(4) 


Ovary 


0.205(1) 


0.209(2) 


0.497(6) 


0.373(4) 


0.435(5) 


0.317(3) 


0.831(8) 


0.680(7) 



Numbers into brackets show stability values sorted out for each tissue from 1 (most stable) to 8 (less stable). The stability values of the most stable gene 
for each group are marked in bold. 



1670 



Zhang etal. (2013) Asian Australas. J.Anim. Sci. 26:1665-1671 



Bestkeeper were slightly different, the two groups of genes 
showing the best (18S, TBP) and the poorest (GAPDH, 
EEF1A2) stability were the same (Table 3). This implies 
that GAPDH and EEF1A2 are not suitable to be used as the 
reference for quantitative analysis of gene expression in 
goat tissues. In addition, the expression of the most 
frequently used reference gene ACTB was not sufficiently 
stable in the present study. The optimal combination of 
reference genes obtained by using geNorm was 18S and 
TBP, but the paired variation coefficient of them was higher 
than 0.15 (V 2 /3 = 0.234), indicating that they were not 
sufficient to be used as the criteria for calibrating the 
quantitative analysis for gene expression in goat tissues, and 
other candidates of reference genes were required (Figure 2). 
The paired variation analysis indicated that the variation 
coefficient showed 0.183 (V 4/5 , the smallest among all 
variation coefficients) for HMBS when its stability values 
were added. Using more than three reference gene 
candidates in combination may increase the level of 
difficulty and introduce errors into the experiment, and may 
also increase the cost of the experiment (Vandesompele et 
al., 2002; Schmid et al., 2003; Lardizabal et al., 2012). We 
therefore recommend that in the present study, the 
combination of three reference genes (18S, TBP, and 
HMBS) should be used for the calibration of the results 
from quantitative analysis of gene expression in goat tissues. 

The present study was carried out to identify suitable 
reference genes (combination) in ten different tissue types 
sourced from Boer goats for the first time by using qRT- 
PCR technique. The results showed that the frequently used 
reference genes in other mammals (pigs, cattle, sheep and 
others) may not be suitable for goats. The expression of 18S, 
TBP, and HMBS was the most stable in the ten goat tissue 
types, and their combination was the most suitable to be 
used as the reference for quantitative analysis on gene 
expression in goat tissues. Our future studies will focus on 
consolidating the best optimal reference genes identified in 
the present study. We will use the results to calibrate the 
quantitative analysis of major functional genes in goats in 
order to obtain accurate, reliable and significant data. The 
data will provide a solid theoretical basis for illustrating 
molecular regulatory mechanisms of economically 
important traits in goats and accelerating variety/breed 
improvements in goats. 

CONCLUSION 

In conclusion, this work evaluated the stability of eight 
genes in different goat tissues showing that they could be 
used as reference genes. Stability values reflect that tissue is 
an important variability factor and it must be taken into 
consideration in the experimental design. It is recommended 
the use of three reference genes: 18S, TBP, and HMBS in 



studies which include multiple tissues. For studies in a 
specific tissue, the most stable reference gene varies 
between the tissues studied. In stomach, small intestine and 
ovary it is recommended ACTB. In heart and spleen the best 
option is 18S, while in uterus and lung we recommend 
using HMBS. If the study is focused in liver, we encourage 
using TBP, but if working with kidney and muscle, the most 
stable gene to be used as reference gene is HPRT1 and 
GAPDH, respectively. 
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